Ultra-discrete Optimal Velocity Model: a Cellular- Automaton Model for Traffic Flow 

and Linear Instability of High-Flux Traffic 
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In this paper, we propose the ultra-discrete optimal velocity model, a cellular-automaton model 
for traffic flow, by applying the ultra-discrete method for the optimal velocity model. The optimal 
velocity model, defined by a differential equation, is one of the most important models; in particular, 
it successfully reproduces the instability of high-flux traffic. It is often pointed out that there is 
a close relation between the optimal velocity model and the mKdV equation, a soliton equation. 
Meanwhile, the ultra-discrete method enables one to reduce soliton equations to cellular automata 
which inherit the solitonic nature, such as an infinite number of conservation laws, and soliton 
solutions. We find that the theory of soliton equations is available for generic differential equations, 
and the simulation results reveal that the model obtained reproduces both absolutely unstable and 
convectively unstable flows as well as the optimal velocity model. 
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I. INTRODUCTION 



Started out as research on vehicular traffic, traffic flow 
has connected with a wide range of social problems such 
as traffic jam [H,[3,lll> evacuation in emergency 0,@, and 
biological systems Also, traffic flow is a good subject 
to apply mathematical developments for 0, [1, H, lioj ] . 
Highway-like traffic, intended in this paper, is modeled 
as a one-dimensional system in which a number of parti- 
cles move in the same direction interacting one another. 
One may regard it as compressible fluid from the macro- 
scopic viewpoint or as a many-body system of driven par- 
ticles from the microscopic viewpoint. Accordingly, there 
appear a lot of models [l], 0, EE E2, EH EH; some are 
coupled differential equations, and others are cellular au- 
tomata (CA). CA models are fully discrete, i.e., not only 
independent variables but also dependent variables take 
integer values, and hence it is quite suitable for computer 
simulation. 

What makes traffic flow a distinct subject from tra- 
ditional physics is that particles in traffic flow have an 
asymmetric interaction. Each particle moves mainly un- 
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der the influence of the particles in front but in contrast 
hardly does under that of the following ones. For this 
idea, in earlier works [TH EH, the car-following models 
were introduced as 

i n (t + T)=F{h n (t)) (h n (t) = x n+1 (t) - xjt)) (1) 

where F is a function, r is a constant and x n (t) de- 
notes the position of the nth particle at time t; ac- 
cordingly, h n (t) means the headway to the particle in 
front, i.e. (n + l)th one. It is remarkable that the car- 
following models introduce two fundamental ideas: one 
is time delay for the particle to respond to the change of 
traffic situation, and the other is a function which pre- 
scribes the optimal velocity for the distance. Namely, 
Eq. JTJ) means that each particle tends/intends to ad- 
just its own velocity to the optimal velocity prescribed 
by F, which requires a delay r. It is noted that the 
car-following model has an exact solution if one chooses 
F(h) = V [l - exp(( 7 o/Vo)(ft - L ))], where ViuL , and 
7o are constants determined empirically [la . Il6j |. An- 
other choice will be seen below. For the latter purpose, 
we refer to the function denoted by F in Eq. ([T]) as the 
optimal velocity (OV) function. 

If t is small enough to make the approximation x n (t + 
t) ps x n (t) + Tx n {t) in Eq. l[Tj). we have the optimal 
velocity (OV) model [l2| 



x n (t) = a[F(h n (t)) -x„(t)}, 



(2) 



where a (= 1/r) means the sensitivity of the particle to 
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traffic situation. The OV model successfully reproduces 
the linear instability of uniform flow in which every parti- 
cle maintains the same distance from the particle in front. 
In particular, high-flux traffic appears as a uniform flow 
with middle headway. (Note that the flux is defined by 
the product of density and velocity.) The stability of the 
OV model changes depending on both the sensitivity and 
the uniform headway, and there exists a critical line at 
which the uniform flow changes from unstable to stable. 
If one chooses 



F(h) = tanh(/i - c) + tanhc (c > 0), 



(3) 



then one finds the critical point (a c ,b c ) — (2, c) (a c and 
b c denoting the critical values of the sensitivity and the 
uniform distance, respectively) at which the flux reaches 
its maximum value [12j. In what follows, we refer to the 
OV model with Eq. ([3]) as the OV model, and the car- 
following model with Eq. ((3J as the delay OV model. It 
is remarkable that the delay OV model also has an exact 
solution expressed by an elliptic function [i| [l?], HH • 

In [l^|, they show that, around the critical point men- 
tioned above, the OV model reduces to the modified 
Korteweg-de Vries (mKdV) equation (20j 



br = b 



xxx 



(b 3 )x 



(4) 



(subscripts X and T denoting partial differentiations in 
them). Actually, in the scaling limit: 



a = 2(1 - e 2 



-2eb(X,T), 
4 (5) 



X = 2e(n-t), T=je% 

where the scaling parameter e tends to zero, one obtains 
Eq. from Eq. @ as the lowest order approxima- 
tion. The mKdV equation is one of well-known soliton 
equations, and it admits kink solutions [2l| . Kink solu- 
tions show a shock wave with some plateaus, and it can 
describe a traffic jam propagating to the upper stream. 
In [191, they also claim that a kink solution is selected, 
independently of the initial condition, by the next-order 
correction and one therefore observes a kink solution at 
the critical point in numerical simulation. 

Soliton equations possess quite rich mathematics such 
as a series of soliton solutions, hierarchy of equations and 
an infinite number of conserved quantities [2(| A 
iV-soliton solution contains N solitary wave packets, i.e. 
solitons, each of which behaves like a particle. One can 
import this solitonic nature directly from soliton equa- 
tions to CA by using the ultra-discrete method 22]. The 
ultra-discrete method is a nonanalytic limiting procedure 
to discretize the dependent variable of difference equa- 
tions. The obtained CAs with solitonic nature are called 
the soliton cellular automata (SCA) or the box and ball 
systems (BBS) [22],[23|]. It is, however, rather complicated 
to obtain a SCA from a soliton equation since one has 
to discretize the soliton equation with respect to both 
time and space variables in the first stage. In the fol- 
lowing discussion, we call equations semi-discrete if the 



space variable is discrete, and full-discrete if both time 
and space variables are discrete. (Note that, each for 
continuous soliton equations, there often exist a number 
of different semi-discrete and full-discrete equations that 
reduce to the same continuous equation.) 

In this paper, we propose the ultra-discrete optimal ve- 
locity (udOV) model derived from the OV model by the 
ultra-discrete method, and exact solution of the model. 
Also, we present numerical simulation results to assess 
the validity of the CA model for traffic flow. The present 
result suggests a basic relation between deterministic 
traffic-flow models and the soliton equations, which con- 
tribute to theoretical studies of traffic flow. 



II. ULTRA-DISCRETE OPTIMAL VELOCITY 
MODEL 

In order to obtain the udOV model, we focus on the fol- 
lowing facts: the OV model reduces to the mKdV equa- 
tion in an appropriate scaling limit 1 191 . and the mKdV 
equation is transformed into a SCA [24|. In addition, it 
is crucial that the delay OV model and the semi-discrete 
mKdV (sdmKdV) equation given in [24| are to meet if 
one supposes a traveling-wave solution. 



A. Ultra-discrete mKdV equation 

First, we shall review the ultra-discrete process of the 
mKdV equation given in 24] , and see that it is natural to 
choose the delay OV model in the ultra-discrete process. 

To begin with, we consider the full-discrete mKdV 
(fdmKdV) equation proposed in [25| (see also pel]): 



1 



i 



„*+! 



1 



(6) 



where S and a are real parameters, and discrete time 
and discrete space variables are denoted respectively by 
super- and sub- scripts, t and j. (Note that we also use 
subscripts for partial differentiations as far as it does not 
cause confusion.) Taking the limit S — ► after substitu- 
tion Vj = rj(—St) in Eq. ([6]) yields the sdmKdV equation 
or the modified Lotka-Volterra (mLV) equation 



rj =r J {l + ar j ){r j+1 



(7) 



where the time variable becomes continuous and we de- 
note by dot differentiation in time. Again, as for a scaling 
parameter e which is introduced as 

1 



2a 



f -les(X,T), 



(8) 



T = —t, 
3 ' 



taking the limit e 

st + 6as 2 sx 



0, we recover the mKdV equation 
1 



4a 



sxxx 



0, 



(9) 



3 



where s — s(X,T) depends on continuous variables in 
both time and space, and subscripts denote partial dif- 
ferentiations in these variables. Thus, one sees that it is 
plausible to adopt Eq. as the fdmKdV equation and 
Eq. ([7]) as the sdmKdV equation. 

The ultra-discrete method means the transformation 
of the dependent variable in full-discrete equations by 
the following formulas: 



.4 






■ exp -j 


e 




.4 


B 


e 


■f exp — 



B, 



(10) 



where A and B are called ultra-discrete variable. Namely, 
the ultra-discrete method is done by taking the limit 
where the dependent variable diverges as the scaling pa- 
rameter e tends to zero. One may formally consider that 
the ultra-discrete method is to replace multiplication with 
•plus, and plus with max as well. It is apparent that the 
commutative, associative, and distributive laws hold. 

Before applying the ultra-discrete method, we intro- 
duce another dependent variable in Eq. © as 



3 1 + avi 



(11) 



Thus, we have another form of the fdmKdV equation 
_ x 1 + (S - a)v%\ 1 + (J - a)v)_ x 



3 1 - oSftl 



1 



av)_ x 



(12) 



Then we shall apply the ultra-discrete method for Eq. 
(TT2")) . The ultra-discrete variables are introduced for both 
the continuous dependent variable and the parameters as 



u.- = exp — , 6 = exp ■ 
J e 



-D 



-A 



exp ■ 



(13) 



where D and A are integers. Taking the limit e — > +0, 
one obtains the ultra-discrete mKdV (udmKdV) equa- 
tion 



Vj +1 + max(0, Vgl - D, Vgl - A) 
- max(0, V#i - A) 
= Vj + max(0, Vf_i - D, Vj_ x - A) 
-max(0,V£_i - A). 



(14) 



As far as we have integers as the initial value, Vj takes 
integer values. Especially if < A < D then Eq. (TTi|) 
reduces to Vj +1 = Vj. 

The udmKdV equation realizes the box and ball sys- 
tem with a carrier (BBSC), which is an extended BBS. 
However, we omit further details since it is not neces- 
sary for our later discussion. Refer to [U, [27[ for more 
information about the BBSC. 



B. The delay Optimal Velocity model 

As described in the Introduction, it has often been 
pointed out that there is a close relation of traffic-flow 
models to the soliton theory [1 [H, [H El . Now, follow- 
ing [l8| we show the direct connection between the delay 
OV model and the mLV equation, which has already ap- 
peared as a semi-discrete mKdV equation in Eq. (JT]). 

We start with the delay OV model in the headway 
representation, i.e., 

h n (t + T)=F(h n+1 {t))-F(h n (t)). (15) 

(See Eqs. (JTJ) and ©■) Changing the variable, 

g n = tanh(/i„ - c), (16) 

we have 

g n (t) = (l-gl{t))(g n+1 (t-r)-g n (t-T)), (17) 

where we transfer r to the right hand side for the lat- 
ter discussion. Then, if one assumes the traveling-wave 
solution 

S»(t)=G(0) (0 = i + 2nr), (18) 
Eq. (JT7J) yields 

G'(<t>) - (1 - G 2 (</>)) (G(0 + t)-G(0-t)). (19) 
Meanwhile, the change of the variable in Eq. ([7]), 



yields 

- 4a ^ = (1 - *j) ~ f j-0- 
Consequently, if one assumes 



rj = Rty) 



(20) 



(21) 



(22) 



Eq. ([21]) coincides with Eq. (fT9|) , 

The traveling-wave solution Eq. (fT8|) means a traffic 
jam propagating to the upper stream at phase velocity 
— l/2r. It is remarkable that this solution is numerically 
observed and also it appears independently of the initial 
condition. That means that the traveling-wave solution 
is stable to perturbations and is manifested as the sta- 
tionary state in the long time limit. Some exact solutions 
of Eq. (|19p have been given by some research groups in 
different ways 0, [13, EH; however, it is recently proven 
that these solutions are all identical [281. 



C. The full-discrete Optimal Velocity model 

Now, we introduce the fdOV model before the udOV 
model, based on the observations above. 
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According to Eq. 1|20[) . we also change the variable in 
Eq. © such as 



2a 



(23) 



This leads to = rj (—St) , and we thus have the fdmKdV 
equation 



4a- 25 



(v 



t+i 



«j)=d (24) 

-(1- ©*+!)(! +5*)^ 



which reduces to Eq. (f2~lj) in i5 — > 0. 

Instead of Eq. (fT8|) . consider the following reduction 
in Eq. (fl~7|) (allowing half integers in the subscript), 



ffn(* - T ) = 9n-l/2{t) 



(25) 



and we can thereby make the equation symmetric in both 
time and space. Then, renumbering the space index (i.e., 
multiplying it by two), Eq. (Ti"7|) reduces to Eq. (f2~Tj) . 
In order to obtain the fdOV model from Eq. (|24p . we 
assume the following reduction as the full-discrete analog 
ofEq. Ip5): 



1-1/2' 



(26) 



where we denote by u l n the dependent variable for the 
fdOV model and m is a positive integer. 

Comparing Eq. (fill) with E q- p2]K one finds 



(27) 



(t = -t'/(4a), n = j/2). 



Then, since —5 is thought to be the unit of time in Eq. 
(|2~4")1 . we introduce the unit of time 7 as 7 = (5/ (4a) in 
Eq. (TTTD, and let 



g„(7i) and 



7-/7. 



(28) 



D. The ultra-discrete optimal velocity model 

We introduce several variables so far; in particular, 
the fdOV model Eq. (|3"0)) is presented as an equation 
for an unphysical variable. In contrast with the case of 
Eq. (|30[) j one may not find a simple equation for the full- 
discrete variable corresponding to h n , i.e., a fdOV model 
corresponding to Eq. (|15p . However, we can obtain the 
udOV model corresponding to Eq. ([15]) as seen below. 
Note that 5/(4a) = t. 

From Eqs. (33]), (H3J), and ([21, one finds 



av) 



1 



1 



V* 
tanh ■ 



.4 



2c 



(31) 



Then, in view of Eqs. (31]), (2S]), and ([23, Eq. (JSU) sug- 
gests that the ultra-discrete variable il* and parameter 
C for the udOV model should be taken as 



gn{rt) — u n — tanh - 



Hi 



c 



2c 



namely, 



h n {rb) 



Hi 

2c 



and c = 



C_ 

2c 



(32) 



(33) 



Accordingly, the ultra-discrete variable H\ has an ex- 
plicit meaning: it is the headway discretized with unit 
length e. 

Here, as well as in Eq. (Tl2l . we introduce as 



exp ■ 



Hi — C 



(34) 



Together with Eq. 



ui +1 



1 + (1 - 4rR +1 



1 + u l 



we thus transform Eq. ([3H)) into 

^ 1+(1 ~ 4 : ) 5' 1 . (35) 



n+l 



Thus, Eq. (p6|) certainly reduces to Eq. (|25|) in the limit 
7^0. Accordingly, we shall consider the following cor- 
respondence in Eq. 



u j+i — "n+1/2 



t-m+l 



W J-1 



7/' — 1l t ~ m 

"n-1/2 ~ u n ' 



(29) 

and thus obtain the full-discrete equation from Eq. (fT7|) . 
Suppose the approximation mentioned in the Introduc- 
tion: t is small enough for the reduction of the delay OV 
model to the OV model. Then we should take m — 1 
which means r = 7. Namely, r is also regarded as the 
minimum unit of time. Consequently, we obtain the full- 
discrete optimal velocity (fdOV) model as 



1 - 2t 



t+i 



<)=(i-<)(i + < +1 K 



+1 

* "\„,t-i 



(30) 



(i-< +1 )(i + <) u : 



In fact, Eq. (|3T)]) is a second-order difference equation 
containing three sequential points of time: t— 1, t, and £+ 
1, and which is consistent with the OV model, a second- 
order differential equation. Also, Eq. (|30[) reduces to Eq. 
if one assumes Eq. (|26p . 




FIG. 1: The increment H^ +1 — H l n is illustrated for each 
region. Regions (A*), (B*), and (C*) divided by dashed 
lines are identical to those in Eq. (|38|l . The superscripts + 
and — show the sign of the increment in each region. The 
increment is continuous, and in particular it takes the value 
of zero on the gray lines. Note that — T < H^ +1 - ff* < T. 
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Suppose < t < 1/4 and 1 — 4r — > +0 when apply- 
ing the ultra-discrete method, we then introduce another 
parameter T as 

r =i(l-exp^) (T>0). (36) 

Consequently, taking the limit e — > +0 after substitution 
of Eqs. (f3"4"| and (|3"6")l into Eq. (|3"B")1 . we obtain ifte ultra- 
discrete optimal velocity (udOV) model 

H' n +1 + max(0, H l n+1 - C - T) - max(0, - C) 

= U l n + max(0, i?*" 1 - C - T) - max(0. H^ 1 - C). 

(37) 

(Note that, as far as T and C are integers, if* also takes 
integer values.) 

We may also express the udOV model as 

—Hi~ l 



- Hi = 



C + T (A + ) 



C-T 



H l n+1 - C 

TTt-l 
— tl,„ 



n n+l 

±T 




f c 

rxt-1 



(A") 
(B+) 

(B-) 
(C ± ) 



(38) 



where regions (A 1 * 1 ), (B ± ), and (C^) are illustrated in 
Fig. Q] This expression explicitly shows how changes 
depending on H^ +1 and H^ 1 . More precisely, each par- 
ticle regards its headway at the previous time step, and 
then decides its motion taking the present headway of 
the front particle into account. 



III. SIMULATION OF THE UDOV MODEL 
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FIG. 3: h = C + T - 1 = 6. Hi at (a)t = 0, (b)t = 10, 
(c)t = 70, and (d)t = 180. (e)time evolution of Ji"^ (n = 40). 



"A 



In this section, we present simulation results of the 
udOV model with the open boundary conditions, i.e., 



FIG. 4: h = C + T = 7. ff* at (a)t = 0, (b)t = 10, (c)t = 60, 
and (d)t = 170. (e)time evolution of Hi (n = 30). 



FIG. 2: h = C = 4. Hi at (a)t = 0, (b)i 
and (d)t = 180. (e)time evolution of ff„ (n = 



= 15, (c)t = 70, 
40). 



there are an infinite number of particles numbered with 
integers. In order to investigate the instability to per- 
turbations, we take several uniform flows each with a 
constant headway (denoted by h below) at t = — 1 and 
then impose a perturbation at t = 0, i.e., the headways 
of several preceding particles are displaced at most by 
±1. We have C = 4 and T = 3 by reference to the real 
data given in (the length of a car)=5m, a = 2s -1 , 
and the inflection point of the OV function (20m, 15m/s). 
In the following simulations, we observe N = 100 parti- 
cles and — 1 < t < 200, which is an enough time for the 
present system to reach a steady state. Then, we let 10 
preceding particles have disturbed headways. 

Figures show Hi (n — 1, . . . , N) for four points 
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of time (t = 0, 10, . . .), and if* (t > —1) for an arbi- 
trary nth particle. Since the OV function is assumed to 
be a monotonically increasing function, a large headway 
means that the particle moves at a high speed. In the 
early stage, we see some features common in the figures: 
the small perturbation initially imposed grow quickly and 
particles then have a headway either smaller than C or 
larger than C + T (to be precise, the headway may take 
these values); moreover, the headways do not take values 
smaller than C — T or larger than C + 2T. In the sub- 
sequent stage, we observe distinct patterns according to 
the uniform headway h as follows. (i)h < C (Fig. [5]): the 
disturbance is amplified and then splits into two oppo- 
site parts, i.e., a free flow region and a traffic jam. The 
free flow region, increasing in size, moves away to the 
upper stream, and eventually all the particles numbered 
n < N are caught in the traffic jam. (ii)C < h < C + T 
(Fig. [3]): a free-flow region, accompanied by a traffic jam, 
moves to the upper stream as well as in (i). Then, in this 
case, an oscillatory pattern follows, in which particles 
stop/decelerate when catching up with the front particle 
and start/accelerate again as the headway becomes large 
enough. Such a state is often referred to as a stop-and-go 
state [H- (iii)ft >C + T (Fig. gj): after the disturbance 
grows into a traffic jam, it moves to the upper stream 
increasing in size. However, particles having got through 
the traffic jam take a larger headway than before, and 
the traffic jam passes out of the frame after all. In [3l| . 
such a uniform flow is referred to as convectively unsta- 
ble in distinction to the absolutely unstable uniform flow, 
which appears in (i). 

The above numerical results can be roughly explained 
by Fig. [T] It is apparent that the system is stable when 
in the regions denoted by zero. Particles therefore tend 
to have headways larger than C + T or smaller than C. 
Moreover, as far as the present initial condition is chosen, 
we conclude that C - T < H l n < C + 2T since —T < 

H n +1 ~ H n< T - If h takes the value from C to C* + T, 
then it always perturbs the steady state with a constant 
headway. 



IV. EXACT SOLUTION OF THE UDOV MODEL 

In this section, we derive an exact solution of the udOV 
model from the one-kink solution of the udmKdV equa- 
tion; in particular, the ultra-discrete one-kink solution 
describes a shock wave traveling with a constant veloc- 
ity. 



A. One-kink solution of the udmKdV equation 

We start with the one-kink solution of Eq. ((9]): 

S (X,T) = ±^tanh(x + i-T), (39) 



which is immediately obtained by assuming a tanh so- 
lution. (The TV-kink solution is also given in (s^l-) In 
the same manner, we obtain the semi-discrete one-kink 
solution of Eq. (j2"Tj) as 



tanhi 
2a 



■ tanh ^ or/ 



tanh < 
2a 



-t 



(40) 



Consider the series expansion: tanh a = a — a 3 /3 + 
in Eq. ([40|) . and we have 



. t anh((i-l)a+±a 3 - 



1 



2n 



(41) 

which suggests Eq. ([5]). 

In order to obtain the full-discrete one-kink solution 
of Eq. ([6]), we use Hirota's method. Hirota's method 
consists of two steps: the first is to transform the soli- 
ton equation into a system of bilinear equations, and the 
second is to solve the bilinear equations. The depen- 
dent variables of the bilinear equations are called the tau 
function. As far as soliton solutions are concerned, the 
tau functions have a simple form and we can hence ob- 
tain soliton solutions. Now, we use the full-discrete tau 
functions given in (33j . i.e., 



0- 



A t+1 

-W 1 



(42) 



where (3 is a constant. Then, following Hirota's method 
we assume the tau functions as 



a) = 1, 



(43) 



where K and L are the parameters to be determined, 
as well as /3, by substituting into Eq. (J6j) . We thereby 
obtain the kink solution which reduces to Eq. (|40[) in the 
same limit as Eq. © reduces to Eq. ([8]): 



^l + L'iP" 1 



v'=f3- 



1 + L l Ki 

where K and L satisfy the dispersion relation, 

(A - l)K 2 L 2 + (K 2 - 2\K + 1)L + A - 1 = 0, (45) 

where A = 6/ a. Calculation to show that Eq. (|4"4f re- 
duces to Eq. (|4"0)) in the limit <5 — i> is relegated to 
Appendix „ 

Before applying the ultra-discrete method, recalling 
Eq. (JTTJ) we have from Eq. (JUJ) 



ILK - 1 K + L t K j 
a K — 1 1 + L t + 1 Ki+ 1 



(46) 



In applying the ultra-discrete method for Eq. (|46|) . we 
assume that v* > is always true. Hence, letting in Eq. 

(USD 



K = exp ■ 



P 



exp • 



Q 



LK — 1 
K — 1 



B 



exp ■ 



(47) 
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where B is also to be determined by K, one obtains the 
ultra-discrete kink solution as 

V? = A+P+max(P, Pj+Qt) -max(0, (j+l)P+(t+l)Q). 

(48) 

In addition, we apply the ultra-discrete method for Eq. 
(l4"5)l and for the definition of B given in Eq. (|4"T)) . and 
thus have the ultra-discrete dispersion relation for Eq. 
(® as 



-(Q-max(0,A-D)) = max(P+Q, 0)-max(P, 0). (49) 



Then, B is equal to each side of Eq. 



B. Exact kink solution of the udOV model 

In the following discussion, we restrict ourselves to 
K > 1 and LK > 1 without loss of generality. (In 
fact, due to the symmetry of Eq. © and that of 
Eq. (0S]), we have v^K' 1 , L~ l ) = Kv){K,L)\ where 
we denote Eq. (|4"4"|) by Vj(K,L) to emphasize the pa- 
rameters therein.) Accordingly, we have P > and 
B = Q = -max(0, A - D). If A < D, then Eq. (fig]) 
becomes V* +l — Vj, i.e., time dependence vanishes. If 
A > D, we have 

V* = A+Q+max(P,Pj+Qt)-max(0, (j+l)P+(i+l)Q), 

(50) 

where Q = D — A. 

Assume in Eq. ([T4")) that A > D and the traveling- wave 
solution, 



v/=nvo &=pj+Qt). 

Then we have 

+ Q) + max(0, V(?/> + P + Q) - £>) 
- max(0, V(tp + P + Q)-A) 
= V{ip) + max(0, V{ip - P) - D) 
- max(0, V{ip - P) - A). 



(51) 



(52) 



Meanwhile, Eq. ([BD]) becomes 

V{ip) = D + max(P, ?/>) - max(0, V + P + Q). (53) 



FIG. 5: An exact solution of the udOV model given in Eq. 
J56} with (7 = 4, T = 3; (a)t = -100, (b)t = -99. This 
shows that traffic jam (region of small headways) propagates 
to the upper stream, taking the shape of (a) at even time, 
and (b) at odd time. 



As well, letting 

Hi = H{4>) ((/) = 2kn + ujt), 

we have from Eq. (f3"T)) 

H(4> + uj) + max(0, H(<f) + 2k) — C — T) 
- max(0, H((f> + 2k) - C) 
= H{4>) + max(0, H(cj) - w) - C - T) 
-max(0,P(</--w)-C). 



(54) 



(55) 



Compare Eq. ([55)) with Eq. ([52^) . and we finally find an 
exact solution of Eq. ([37)) which presents a shock wave: 



H(</>) = C + T + max(T, cj>) - max(0, <j> + 2T), (56) 

where = (2n + t)T and C, T > 0. Figure O shows a 
kink solution of Eq. (|3T[) given in Eq. ([56]) . 



V. SUMMARY AND CONCLUDING REMARKS 

In this paper, we propose the ultra-discrete optimal 
velocity model, a cellular-automaton model which inher- 
its the properties of the OV model such as the linear 
instability, and moreover have an exact solution of the 
model which describes a traffic jam propagating to the 
upper stream like a shock wave. Since the OV model 
is defined by a nonlinear differential equation, ordinary 
discrete schemes do not work to obtain the correspond- 
ing difference equation. Meanwhile, it is well-known that 
the OV model, around the critical point, reduces to the 
mKdV equation, a soliton equation, in a scaling limit. 
Therefore, we apply the theory of soliton equations to 
the OV model so that the CA model obtained maintains 
the relation between the OV model and the mKdV equa- 
tion described below. 

Figure [S] shows our scheme of the present study. At 
the beginning, we pay attention to the fact that if one 
considers traveling-wave solutions, the delay OV model 
Eq. (fTT)) coincides with the sdmKdV equation or the 



mKdV 



OV 



continuous 



(9)- 



(5) 



"(2) 



semi-discrete 



(8) 

(7) = (21). 

(20) 



-(17) = (15)- 
(16) 



"(I) 



(13) 



ultra-discrete (14) 




(37) 



FIG. 6: Scheme of the present study. The numbers indicate 
the equations in previous sections. 



8 



modified Lotka-Volterra equation Eq. ([7]) , as well as that 
the OV model Eq. ^ reduces to the mKdV equation 
Eq. (|9|) in a scaling limit Eq. (|5|). We therefore adopt 
the delay OV model Eq. (jTTJ) as the semi-discrete OV 
model and then have the fdOV model Eq. (|30|) from the 
fdmKdV equation Eq. |6]); we need to establish the fdOV 
model before the udOV model. Thus, we have the udOV 
model Eq. ([57)1 and an exact shock solution, and one sees 
the validity of the soliton theory for study of traffic flow. 

Simulation results show that the udOV model has a 
transition region where it is difficult for a uniform flow 
to maintain the headway under perturbations. (It is 
also remarkable that a deterministic traffic-flow model 
has such instability to perturbations.) Moreover, we find 
that the system have three distinct states: a traffic jam 
(absolutely unstable), an oscillatory pattern (stop-and- 
go), and a free flow (convectively unstable) according to 
the uniform headway. 

In the present paper, we mainly focus on the theo- 
retical study and do not have the necessity to have the 
model in position representation; however in some cases 
the spatio-temporal patterns are also required. It will be 
discussed in subsequent publications. 



sity of Tokyo. 



APPENDIX: KINK SOLUTION OF THE 
FDMKDV EQUATION 

We show that the full-discrete one-kink solution Eq. 
(|44|) surely reduces to the semi-discrete one-kink solution 
Eq. (gDJ in the limit S — > 0. First, we solve Eq. (gSJ in 
L and take one of the two solutions such that (3 in Eq. 
(144p does not diverge as 5 — > 0: 



L 



K 2 - 2XK + 1 + (K - l)y/(K + 1) 2 - AXK 



-1 



K -IS 



K + 1 a 
Consequently, we have 



2(1 -X)K 2 
0(S 2 ). 



L i/s 



K - 1 



exp ■ 



a(K+l) 
Next, we transform Eq. (|4"4"| into 



— 0). 



(A.l) 



(A.2) 
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1 (LK- 1)(K+ 1) 



2a 



LK 2 — 1 



1 



K - 1 L l K' J - 1 



Let K = e 2a , and we finally have 
-t/s 1 



K + 1 L*KJ + 1 



<5^o 2a 



tanh a ( tanh a 
tanh aj 1 



2a 



2a 



(A.3) 



(A.4) 
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